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ABSTRACT 

We present multi-epoch medium-resolution observations of two M4.5 candidate members in 
the halo of the ^8 Myr rj Chamaeleontis open cluster. Over six months of observations both 
stars exhibited variations in their Ha line profiles on timescales of days to months, with at 
least one episode of substantial activity attributable to accretion from a circumstellar disk. We 
derive an accretion rate ~10 -87 M© yr -1 for this event, with a rate of ^10 -10 ' 6 M© yr _1 in 
quiescence. Episodic accretion like that observed here means existing surveys of accreting 
Weak-lined T-Tauri Stars in young clusters are likely incomplete and that gas dissipation 
timescales calculated from the fraction of accreting objects are underestimates. 
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■ 1 INTRODUCTION 

t-H ■ The open cluster rj Chamaeleontis is one of the closest (94.3 pc) 
and youngest (~8 Myr) stellar aggregates in the Solar neighbour- 
. hood. Nearby, isolated groups such as rj Cha are ideal laborato- 
' ries for investigating the dynamical evolution of young star clus- 
ters, in particular the influence that dynamics have had on the evo- 
lution of protoplanetary disks. Young clusters show a steady de- 
cline in the n umber of stars having disks and signatures of a ccre- 
tion with age JMohantv et a l. 2005; Ja vawardhana et al .1120060 . By 
an age of ~5 Myr, 90-95% of all young cluster members have 
stopped accreting material at a significant rate, yet ~20% of ob- 
jects retain enough dust in their disks to produce a mid-IR excess. 
The mechanis m responsible for these two different timescales is 
still uncertain (Fedele et al. 2010). rj Cha has already been shown 
to have both an excess of stars with accretion and st ars with de- 
tectable disks compared to clus ters of similar age dHaisch et al.l 
l200ll : ISicilia-Aguilar et alj|2009l) . The cluster also appe ars to have 
an In it ial Mass Function d eficient in low-mass objects iLvo et all 
l2004h . iFedele et alj (2010) is the latest study to point out that dy- 
namical evolution has probably dispersed a large fraction of the 
low-mass members to radii beyond that currently surveyed (see also 
Moraux ct al Murp hy et al ■l2010h . Depending on the proper- 

ties of the dispersed members, it is possible that the current mem- 
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bership may be biased towards IR-excess and strong Ha emitting 
sources. By investigating the disk and accretion properties of any 
new dispersed members of rj Cha we can hope to gain a more un- 
biased view of the cluster as a whole, as well as addressing any 
influen ce dynamical interact ions have had on disk evolution. 

In iMurphv etafl boil hereafter Paper I) we presented the 
results of our search for the putative halo of low-mass objects sur- 
rounding rj Cha. From photometry, proper motions and multi-epoch 
spectra we identified 4 probable and 3 possible members, includ- 
ing two M4.5 stars which exhibited large variations in their Ha 
emission line strengths. These stars, 2MASS 08014860-8058052 
(hereafter 2MASS 0801) and 2MASS 08202975-8003259 (here- 
after 2MASS 0820) are the subject of this letter. 



2 MULTI-EPOCH SPECTROSCOPY 

As part of our continued investigation into rj Cha we obtained 
multi-epoch, medium resolution spectroscopy of 2MASS 0801 
and 2MASS 0820 with the WiFeS instrument at the ANU 2.3 m 
telescope (see Paper I for more details). We obtained 9 observa- 
tions of 2MASS 0801 over 2010 January- June and 13 epochs of 
2MASS 0820. In addition to Ha emission both stars recurrently 
showed He 1 5876 A, 6678 A and Na I D emission, often associated 
with accretion. Details of these observations are given in Table[T] 



2.1 Accretion diagnostics 

To trace any accretion across our observations we analysed the 
time evolution of the Ha equivalent width and velocity width of 
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Table 1. Multi-epoch WiFeS observations. 
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the Ha line at 10% intensity (i>io). Recently iNguven et al.1 (2009) 
have shown that due to the strong dependence of vio on the shape 
of the Ha line profile, the Ca II 8662 A line flux is perhaps a 
more reliable quantitative diagnostic of accretion rate. The Ha uio 
width has a well-established recor d as an accret i on in dicator and 
despite the caveats put forward bv lNguven et al .1 d2009l) . we use it 
here for comparison to previous studies and the fact that our wave- 
length coverage does not extend past 7100 A. Figure Q] shows an 
EW-Dio diagram for our observati ons of 2MASS 0801 and 820, 
with r? Cha member s observed by Javawardhana et aD d2006t) and 
lLawson et all d2004h for comparison. Measurements of the EW and 
vio for each star are listed in Table [T] We estimate an uncertainty 
of ±1 A for the EWs and ±10 km s" 1 in vio, primarily due to un- 
certainties in defining the pseudo-continuum around the broad Ha 
lines at this spectral resolution. Both stars show substantial vari- 
ation in the EW-f io space, with the scatter more pronounced in 
2MA SS 0820. 

lLawson et al. I J2004I) observed the cluster at high-resolution 
and found a similar EW variation of 2-3 x for RE CX 5, 9 and 1 1 
when c ompared to the discovery m easurements of Mamaie k et all 
dl999h . |javawardhana etal] d2006h obtained 3-8 epochs of Magel- 
l&nlMIKE echelle data for 1 1 77 Cha members during 2004 Decem- 
ber - 2005 March. They report variations in EW and «io in some 
stars similar to those seen in 2MASS 0801 and 0820 with 35% mo 
and 50% EW variations in the maximal cases. In Figure[T]we com- 
pare the measurements of these two groups. In general there is good 
agreement between the two sets of measurements, with RECX 9, 
RECX 11 and ECHA J0843.3 all clearly accretors in both stud- 
ies. RECX 7 is a known non-accreting SB2 spectroscopic binary 
dLvo et al.ll2003h. which ex plains the br oad line width but low EW 
in the Javawardhana et al. observations. lLawson et al. I d2004h clas- 
sify RECX 5 as accreting with t?io > 300 km s~ x from their 2002 
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Figure 1. Ho Equivalent Width versus the Ha velocity width at 10% peak 
intensity (vi ) for 2MASS 0820 and 2MASS 0801. Horizontal lines denote 
the minimum EW for CTT stars in the indicated range of spectral types. The 
vertical line at 270 km s separates accreting and non-accreting objects 
IWhite & Basri 2003). Shaded region s show the standard deviation of the 
multi-epoch measurem ents of[ja vawardhana et al. 12006). Comparisons to 
the single epoch data o flLawsonetalJd2004l) are also shown. Lines connect 
the two sets of measurements. For the 2MASS 0820 measurements lines 
connect observations taken on the same night. 



observation. Ijavawardhana et all however failed to detect any ac- 
cretion, with only one of their five observations showing a broad- 
ened Ha line. The blue bump seen in the 2004 December spec- 
tru m could mean RECX 5 is sporadically accreting at low levels 
and lLawson et al.1 observed a strong outburst of 10 -10 Mq yr _1 
accretion in 2002, or the star is chromospherically active and un- 
derwent a period of strong activity during their observations. Chro- 
mospheric activity as a possible mechanism for the variations seen 
in our spectra is discussed in detail in Q2.2\ 

Several sets of 2MASS 0820 spectra were taken during the 
same nights. They reveal a scatter of up to 25 km s in vio and 
several Angstroms in EW, i.e. larger than the instrumental errors. 
While both stars generally lie on the non-accreting side of the 
270 km s _1 criterion defined bv lWhite & Basril d2003l) , they each 
make several excursions into the accreting region of the diagram, 
meeting both the EW and ?;io criteria for accreting Classical T- 
Tauri stars (CTTS). This is similar to the behaviour of RECX 5 
during 2002-2005. The timescale of these excursions appears to be 
of the order hours to days. In the case of 2MASS 0801 both EW 
and vio increase dramatically on Apr 30, before returning to qui- 
escent levels the next night. A similar level of activity is seen on 
2010 Feb 19. Figure [2] shows the evolution of the 2MASS 0801 Ha 
profile during the 2010 April event. We obtained spectra before, 
during and after the event. The pre- and post- event line profiles are 
remarkably similar, with broad wings developing during the Apr 30 
outburst, giving rise to a velocity width of Vio > 320 km s . The 
EW similarly quadruples to —27 A. The central velocity of the line 
profile also evolves with time (see Figure|4j- 2MASS 0820 appears 
to have 3 tiers of activity: a base level at EW ~ —16 A, a higher 
level at —20 > EW > -27 A, then increasing to —40 A and 
Vio = 425 km s" 1 on 2010 Feb 19. 

In addition to WiFeS spectroscopy we obtained a Magel- 
lan/MIKE spectrum of 2MASS 0820 on 2010 May 11. The Ha ve- 
locity profile is shown in Figure [3] with a contemporaneous WiFeS 
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Figure 2. Evolution of the 2MASS 0801 Ha velocity profile during the 
2010 Apr outburst. All profiles have been shifted to a radial velocity of 
zero. The Apr 30 profile has been veiling-corrected as described in the text. 
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Figure 3. MIKE Ha velocity profile of 2MASS 0820 compared to a con- 
temporaneous WiFeS observation. The dotted line shows the R ss 25000 
MIKE spectrum smoothed to the R Rd 7000 resolution of WiFeS. 



iMontes et all Jl998h has shown that in some WTTS the Ha line 
profile cannot be fitted by a single gaussian and two components are 
necessary: a narrow gaussian of FWHM < 100 km s _1 and a much 
broader component with FWHM 130-470 km s -1 , sometimes off- 
set in wavelength from the narrow component. They attribute these 
line profiles to micro-flaring occuring in the chromospheres of the 
stars. Micro-flares are frequent, short duration events and have 
large-scale motions that could explain the br oad wings observed in 
the lines and the residual spectra in Figure [4] Stauffer et al.l J 19971) 
observe a similar effect in several of their ~30 Myr IC 2391/2602 
and 100 Myr Pleiades targets. Both our stars show He I 6678 A in 
emission at their peak Ha levels. Strong He I 6678 A emission is 
an accretion diagnostic as it is generally only pre sent in low-level s 
(<Cl A) in older chromospherically active stars dGizis et al]|2002h . 
While we do detect strong (1.5 A) emission in the April 30 outburst 
spectrum of 2MASS 0801, at all other epochs where we dete ct the 
line it is weak (~0.5 A). Furthermore. iMartm & Ardilal ( 1200 II) de- 
tect strong (1^1 A) He I 6678 A emission during a flare of the old 
M9 dwarf LHS 2065. 

Given the weak He I line strengths generally observed in our 
stars and the simple gaussian-like profiles of the residual spectra 
we do not have strong evidence for ongoing accretion. Chromo- 
spheric activity is a much more likely explanation for the observed 
line profiles and EW/^io variations. Only the February 19 spectrum 
of 2MASS 0820 shows a broad, asymmetric residual characteristic 
of accretion. Multiple components are present at velocities up to 
±300 km s _1 , presumably tracing the ballistic infall of material 
from the inner edge of the disk onto the stellar surface. 



spectrum. The EW and Via of the higher resolution spectrum agree 
well with the WiFeS values (see Figure [T). The velocity profile 
shows some self-absorption in the line center as well as a slight ex- 
cess of blueshifted emission at v ~ —100 km s _1 . To test how our 
lower resolution WiFeS spectra affect measurements of EW and uio 
we smoothed the MIKE spectrum to R ~ 7000 and re-measured 
vio and the EW. In both parameters the smoothed value was similar 
to that measured from the contemporaneous WiFeS spectrum, i.e. 
smaller than the original value. We are thus likely underestimating 
the EW and vio from our WiFeS observations. 



2.2 Ha velocity profile variations 

As well as the bulk quantities of Ha EW and Vio we also investi- 
gated the variation in the shape of the Ha velocity profile over our 
observations. The top panels in Figure [4] show the average quies- 
cent spectrum for each star, constructed using several (7-9) non- 
outburst spectra. The standard deviation of these spectra is also 
shown. 2MASS 0820 shows increased scatter around the mean pro- 
file at v ~ +100 km s _1 (and to a lesser extent at —100 km s _1 ). 
This is due to a variable-strength red shifted component at this 
velocity present in some of the observations. The bottom panel 
shows the variation in the residual spectra - the difference between 
the individual spectra and the mean spectrum. Immediately appar- 
ent are the broad residual profiles, tracing velocities up to ±200- 
300 km s" 1 . The Feb 19 epoch of 2MASS 0820 shows a residual 
velocity profile reaching to ±300 km s - , with four distinct com- 
ponents visible and a large red asymmetry. Velocity shifts in the 
peaks of the residual spectra are present at up to several tens of 
km s . 

Instead of invoking accretion, can chromospheric activity 
also explain the strong and varying Ha profiles? The work of 



2.3 Accretion rates 

iNatta et alj d2004h have shown that a good correlation exists be- 
tween the Ha vio velocity width and the mass accretion rate across 
a large range of masses, 0.04-0.8 Mq. If we assume the Ha emis- 
sion is a result of accretion we can use the lNatta et al.l ( 120041) rela- 
tion to derive the mass accretion rate (M) without the need to per- 
form detailed model fits to the line profile. We derive a quiescent 
accretion rate for 2MASS 080 1 of log M = - 1 1 . 3 ± . 3 M© yr" 1 , 
where the uncertainty reflects the variation in Vio and uncertainties 
in the relation parameters. For 2MASS 0820 the quiescent scatter 
is larger and log M = —10.6 ± 0.5 Mq yr _1 . For the event of 
February 19 we calculate an accretion rate of log M = — 8.7 ± 
0.5 M yr" 1 . This is similar to that derived for ECHA J0 843.3 
from detailed Ha profile modeling by lLawson et al.l d2004l) . The 
INatta et alj d2004l) relation gives an accretion rate 1.3 dex la rger 
for the co rresponding v i o velo city given by lLawson et alj Ac- 
cording to Ngu yen et alj 12009) the large scatter in this average 
relation probably reflects object-to-object variation (possibly due 
to evolutionary effects) rather than the effects of variability on 
the usually single-epoch obse rvations as previously thought (e.g. 
IScholz &Ja vawardhana 2006j). As previously mentioned, the ef- 
fect of our lower spectral resolution means we have likely un- 
derestimated Vio so our derived accretion rates are probably even 
higher. 2MASS 0820 shows a 10-80x jump in ac cretion rate be- 
tween quiescence and its active phase. Using the INatta et alj re- 
lation on th e high-resolution RECX 5 velocities of Lawson et al.1 
d2004l) and Ijavawardh ana et al. I d2006l) results in a similar 20 x 
change in accretion rate. Assuming the scatter in quiescent v%o ve- 
locities is a result of the accretion rate varying with time, the vari- 
ations are similar in magnitu de to those seen in the 10 Myr CTTS 
TW Hva dEisner et alj|20icb a nd lower-mass youn g brown dwarfs 
dScholz & Javawardhani2006l : IStelzer et alj2007l) . 
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Figure 4. WiFeS Ha velocity profiles for the 9 observations of 2MASS 0801 (left panel) and 13 observations of 2MASS 0820. In each plot the top panel 
shows the average quiescent spectrum spectrum and the standard deviation of quiescent spectra around the mean (shaded region). The bottom panel shows the 
variation around the mean for all epochs. Note the broad, multicomponent residual on Feb 19 for 2MASS 0820. 



2.4 Continuum veiling in 2MASS 0801 8058 

The April 30 spectrum of 2MASS 0801 exhibits a Li I 6708 A 
equivalent width smaller than in quiescence. This is the effect of 
continuum veiling, where enhanced continuum emission fills in ab- 
sorption lines and reduces the measured EW. Magnetospheric in- 
fall models predict the presence of veiling from accretion shocks, 
as matter from the ac cretion flow reaches the stellar surface (e.g. 
iMuzerolle etai]|2003h . However, this is probably not the case for 
2MASS 0801. Figure|5]shows the full WiFeS spectrum for April 30 
compared to the average 2MASS 0801 quiescent spectrum. Veil- 
ing is readily apparent as an enhancement of the continuum in the 
6200 A CaH band and a corresponding depression of the contin- 
uum around Ha (due to the normalisation at 6100 A). Using the 
quiescent spectrum as a reference we find the veiling at 6100 A is 
~0.22. The observed flux F was veiling-corrected using the fol- 
lowing prescription: 

Fcorr = F — 0.22 x F 61QoA x B(X, T)/B(6100 A, T) (1) 

where B(X, T) is the Planck blackbody function and T = 2800 K. 
With the correction applied the continuum and Li I 6708 A EW 
once again match the quiescent levels and the only excess is now 
due to line emission from Ha, He I and Na D. The Ha EW is in- 
creased from —20 A to the —27 A given in TableQ] Given the low 
temperature of the veiling blackbody the origin of the veiling emis- 
sion is most likely warm circumstellar dust rather than an accretion 
shock, which generally emits at much higher temperatures (5000- 
20000 K). 2MASS 0801 appears to have an excess of ~0.15 mag 
in J - H and -0.05 mag in H - K s relative to 2MASS 0820 
and other mid-M members which could indicate the presence of 
circumstellar material. The star inhabits a region of colour space 
close to two strong accretors known to possess disks - RECX 1 1 
and ECHA J0843.3 have Class I I SEDs, with significant e xcess at 
wavelengths greater than 3 dSicilia-Aguilar et alj|2009h . 

Reddeni ng (i.e. interstellar dus t) could also explain the colour 
excesses. The Schleg el et al]dl998l) reddening along a line of sight 
to 2MASS 0801 is E(B - V) = 0.4 mag and the star sits on a 
ridge of prominent dust emission in the IRAS maps. It is possible 
that the stellar colours are affected by reddening if the dust lies in 



front of the star. From optical polarisation studies the distance to the 
ridge has been estimated to be ~115 pc dClearv et alJll979T) . 40% 
larger than the 83 pc estimated dynamical distance to 2MASS 0801 
from Paper I. An improved distance to the star will be necessary to 
resolve the reddening issue. 

2.5 Radial velocity variations 

Our multi-epoch observations provide insight into any radial ve- 
locity variations in the two stars. The measurement of velocities 
from the WiFeS spectra is described in Paper I and the resultant RV 
values are given in Table Q] From 260 observations of M spectral 
type velocity standards we have established that WiFeS is capable 
of an RMS RV precision of 2 km s~\ The standard deviation of 
the 13 RV measurements of 2MASS 0820 is 1.9 km s~\ consistent 
with a flat velocity curve. The RV from the MIKE measurement, 
17.4 ±0.5 km s _1 , is also consistent with the mean WiFeS value 
of 18.4 km s _1 within the errors. In contrast, 2MASS 0801 has 
a standard deviation of 3.3 km s _1 from 9 measurements and an 
RMS variation nearly five times larger than 2MASS 0820. 

In Paper I we suggested that in conjunction with its elevated 
position in the cluster CMD the RV variation of 2MASS 0801 sug- 
gested binarity. Any such companion would need to be significantly 
cooler and fainter as we do not detect a second set of lines at any 
epoch. The timescale of the RV variation suggests a period of or- 
der days to weeks with a velocity amplitude of less than a few 
tens of km s _1 . Without high-resolution imaging we cannot con- 
strain the orbital semi-major axis other than to say any companion 
is not resolved by contemporary imaging surveys such as DSS2 and 
2MAS S. High resolution spectrosc opy and adaptive optics imaging 
(e.g.[K6hler & P etr-Gotzensl2002h will be necessary to resolve any 
companion. 

Another possibility is that the RV variation is induced by activ- 
ity or surface features co-rotating with the star. lMartm et alj ( 2006) 
detect a 3.5 km s -1 amplitude, 3.7 hr periodicity in the optical RV 
data of the M9 brown dwarf LP 944—20, which they attribute to ro- 
tationally modulated inhomogeneous surface features. Using a toy 
model with a single spot 100—200 K cooler than t he p hotosphere 
covering ~10% of the stellar surface. lReiners et al.N201oT) can gen- 
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Figure 5. The effects of continuum veiling on the 2010 Apr 30 spectrum of 2MASS 0801. All spectra smoothed by a 10 pixel Gaussian kernel and are 
normalised over the 6050-6150 A region (heavy line). The red dotted line shows the final veiling-corrected spectrum, as calculated by applying Equation[T] 



erate a velocity amplitude of a few km s 1 for stars of 2800- 
5700 K, similar to the signal seen in LP 944-20 and 2MASS 0801 . 
We do not find any correlation between the RV and either the Ha 
EW or v\o so any velocity variation is probably driven by rotation 
not chromospheric flaring. The rotational period of 2MASS 0801 
has not been measured but is likely on the order of days, s imilar to 
the other late-type members of 77 Cha (Law son et alj|200"ll) and the 
timescale of observed RV variation. 



3 DISCUSSION 

Our results show that Ha variability in ~8 Myr PMS stars can 
be substantial on both short (hours-days) and long (months) 
timescales. This variation is probably driven primarily by chromo- 
spheric activity, which can generate broad Ha profiles mimicking 
accretion over short timescales. However, we also have evidence 
for at least one accretion event in 2MASS 0820 which requires 
confirmatory follow-up observations. Additional mid-IR observa- 
tions will be necessary to detect the presence of any circumstel- 
lar disk around the star feeding the accretion. Assuming the duty- 
cycle of episodic accretion is low, single-epoch surveys of accret- 
ing objects, especially in the critical age range 5-10 Myr when in- 
ner disks are being cleared and giant planet formation takes place, 
are likely missing a large fraction of accreting objects. Gas deple- 
tion timescales derived from the fraction of accretors are therefore 
likely underestimates. This may provide another explanation for the 
difference in characteristic timescales between mass accretion and 
dust dissipation found by [Fedele et alj ( bold , 2.3 Myr and 3 My 
respectively), without the need to invoke planet formation and/or 
migration in the inner disk as a possible mechanism for halting 
accretion. A larger survey of the disk and accretion properties of 
outlying 77 Cha members, combined with more detailed investiga- 
tion of the true accreting fraction of PMS clusters from multi-epoch 
surveys is needed to resolve this discrepancy. 
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